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Abstract 

The localization (confinement) of coupled plasmon modes, named as gradons, has been studied 
in metal nanoparticle chains immersed in a graded dielectric host. We exploited the time evolution 
of various initial wavepackets formed by the linear combination of the coupled modes. We found an 
important interplay between the localization of plasmonic gradons and the oscillation in such graded 
plasmonic chains. Unlike in optical superlattices, gradient cannot always lead to Bloch oscillations, 
which can only occur for wavepackets consisting of particular types of gradons. Moreover, the 
wavepackets will undergo different forms of oscillations. The correspondence can be applied to 
design a variety of optical devices by steering among various oscillations. 
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I. INTRODUCTION 



Clusters of metal particles at the nanometer scale have attracted much interest because of 
their plausible applications in nano-optics. 1 The main objective is to overcome the diffraction 
limit, achieving subwavelength confinement of waves by surface plasmons. 2 In the last decade, 
enormous developments have been witnessed in plasmonic waveguides to build all-optical 
circuits. 2,3 Regarding controlled localization of optical signals and guided electromagnetic 
energy at subwavelength scale, there have been various progresses. 4,5 ' 6 Furthermore, Fabry- 
Perot cavity was proposed to confine surface plasmons. 7 Light will be localized when the 
frequency matches the surface plasmon eigenfrequency of the system, while for a wavepacket 
consists of a variety of frequencies, oscillations will occur. 

The most intriguing oscillation is the Bloch oscillation (BO), which is the oscillatory 
motion of a particle in a periodic potential when a constant force is acting on it. It has 
been generally understood that Bragg reflections together with a linear gradient of the po- 
tential can cause Bloch oscillation of a wave of distinct nature (electronic, optical, acoustic, 
or matter wave). 8 The gradient in potential can be introduced by an external field or a 
spatial perturbation of any nature (electric, magnetic, acceleration, or gravitation field). 8 
For example, in optical structures, the optical equivalent of an external field in electron- 
ics can be a gradient in the refractive index 9 or a geometrical variation in waveguides. 10,11 
Many investigations have been devoted to photonic BO in helical (deformed) waveguides, 10 
curved waveguides, 12 waveguide lattices, 13 chirped metamaterials, 14 and other photonic het- 
erostructures. 15 Recently, plasmonic Bloch oscillations (PBO) have been observed in metal 
heterowaveguide superlattices, 16 metal waveguide array structures, 17 and chirped metal- 
dielectric structures. 18 However, the mechanisms of BO, especially the relation between the 
occurrence of BO and the condition of gradient, have not been clearly explained. In this 
work, we aim to explore the mechanisms of dynamic oscillations (including BO) in graded 
plasmonic chains (GPCs) based on our understanding of localized plamonic modes (gradons) 
in previous work. 19,20 

In GPCs, the gradient can be realized by gradually changing the dielectric constant of the 
host, 19,20 the distance between the metal particles, and/or the particle size along a particular 
spatial dimension. 21,22 There exist peculiar confined eigenmodes (or quasi- normal modes) 
which are called gradons in the graded lattices. 19,20 Gradons with different frequencies are 



2 



localized at different position and have different spatial extent in the system. Then one 
would ask the following questions: Can BO occurs in a GPC? What is the mechanism of BO 
in graded optical or plasmonic lattices? Does a gradient guarantee the occurrence of BO? 
Motivated by these questions, we aim to explore the relation between gradon localization and 
the occurrence of BO in GPCs. For this purpose, we shall use a phase diagram to distinguish 
between various plasmonic gradon modes in GPCs. Then through a dynamic simulation, we 
show the evolutions of different "wavepackets" in the position-time domain. We find that 
there exists a clear and interesting correspondence between the gradon localizations and 
the occurrence of BO, breathing- wave-like oscillations, or other kinds of oscillations. We 
explicitly show that only initial "wavepackets" formed by a combination of certain type of 
gradons can undergo BO or breathing-wave-like oscillation. The existence of such kind of 
gradons requires rather strong gradient strength and they are referred as blue-red gradons 
which are well confined in the "internal" region of the system. Conversely, BO and breathing- 
wave-like oscillation cannot occur in GPCs bearing no such blue-red gradons, irrespective 
of whatever initial "wavepacket" . Despite of observing BO in GPCs, we further explore its 
mechanism and find that gradient can not always lead to BO in GPCs. Our results provide 
useful guidelines in designing optical devices of novel functionalities. 



II. MODEL AND FORMULA 



We consider a chain of N spherical silver nanoparticles that is immersed in a dielectric 
host. 19 ' 20 The dielectric function of the metal nanoparticles has the generalized Drude form 23 
e m {uj) = ^ — r]ujp/uj(uj + iT), while the host has a graded dielectric constant varying along the 
chain axis, i.e., x, as €h(x n ) = ei + cx n /l (n = 1,2, N), where x n = ndo is the position 
of the nth nanoparticle, do the interparticle spacing, I = Ndo the total length of the chain, 
and c is the coefficient of dielectric gradient. All the particles are identical in size with radii 
a. The nth and the rath particles are at a center-to-center distance d nm = do\n — ra|. The 
coupled equation for point dipoles reads 20 

G(w) |0) = 0, (1) 

where \<f>) is the column vector of the dipole moments oscillating with frequency u, 
G nn (u) = l//3 n (u) and G nm (u) = -f nm (ra ^ n). Here /3 n (u) = a 3 e h (x n )[e m (uj) - 
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^h(x n )}/[e m (uj) + 2eh{x n )] represents the polarizability of the nth particle, and T nm is the 
dipole electromagnetic interaction between the nth particle and the mth particles, 

f nm = —^(VnVn + A; 2 I) , (2) 

€h\X n ) \Xn X m \ 

where k is the wave vector in the host. We can take the point-dipole approximation in the 
calculation, because the condition a/d < 1/3 is always satisfied. The validity of this approx- 
imation has been proved theoretically and experimentally. 21 To calculate the eigenmodes of 
the system, we assume the external field is zero and take the wave vector k = 0. In the 
quasistatic limit, T nm pa K/ef l (x n )d 3 im , where k = 2(— 1) for the longitudinal (transverse) 
case. 

Actually, the generalized Drude model can be extended from the standard Drude model by 
replacing uj 2 with rjuj 2 /^ and with e^/^, respectively. However, it is hard to determine the 
eigenvalues and eigenvectors by directly diagonalizing G(u), because Eq. (1) is a nonlinear 
complex transcendental equation with respect to the frequency uj. This prevents us from 
constructing a phase diagram which is crucial to understand the connection between various 
gradon modes and the dynamic oscillations. In order to circumvent the difficulties, we 
linearize Eq. (1) with respect to the frequency by using the quasistatic approximation of 
Eq. (2). In this way, we can use standard diagonalization to obtain both the eigenvalues 
and eigenvectors. For noble metal, e.g., silver, T pa 0.005co> p according to the tabulated data 
in the reference. 24 Here we can neglect the loss, i.e., T = in the metal nanoparticles. Thus 
Eq. (1) can be linearized with respect to uj 2 

(F - u 2 I) |0) = , (3) 

where I is an identity matrix and F nm = u 2 — G nm (u 2 ) / G' nn (oj 2 ) . Here the prime indicates 
derivative with respect to uj 2 and u n = 0J p \Jrj/{^ + 2e/j(x n )) is the dipole Mie resonant fre- 
quency of the nth particle at that particular position where the host has dielectric constant 
£h( x n)- Equation (3) represents a very good approximation of Eq. (1) when u pa uo n and it can 
be mapped onto an equivalent elastic chain with mass M n — [£ + 2eh(x n )] 2 /3Xr]u 2 €h(x n ), 
and additional on-site potential U n = M n u 2 — 2K , where K = (a/ do) 3 is a force con- 
stant between adjacent masses. 19,20 Two characteristic frequencies u> c i(n) = \JU n /M n and 
^c2(n) = \J (U n + AK )/M n were defined to construct the phase diagram. 20 The parameters 
used in this work are listed below: N = 100, c = 1.0, e\ = 3.0, £ = 5.45, and rj = 0.73. Thus 
the host dielectric constant varies from e h {l) = 3.0 (left end) to e^(N) = 4.0 (right end). 
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III. PLASMONIC GRADON MODES IN GRADED PLASMONIC CHAIN 



Here we briefly summarize the various localized plasmonic gradon modes and their tran- 
sitions, 20 which are very related to the dynamic oscillations of light wavepackets in the 
graded plasmonic chains. The effective approach of identifying various gradon modes is 
the phase diagram, 25 as shown in Fig. 1, which contains four regions representing four dif- 
ferent kinds of gradon modes. The different phase regions are interfaced by four curves 
denoting, respectively, ^ c i(l) (dashed line), u c i(N) (solid line), cu C 2(l) (dotted-dashed line), 
and u C 2(N) (dotted line) as functions of distance d . In each region, there exists a kind of 
plasmonic gradon modes, which is localized at different position of the system. The mode 
patterns of some plasmonic gradon modes, i.e. square moduli of induce dipoles versus posi- 
tion, are shown in the insets. There exists a critical value do = d c , where u; c i(l) = ou C 2(N). 
For interparticle spacing d < d c , the unbound (extended) modes (U) cover a frequency 
range co> cl (l) < uj < uj c2 (N) (vertical line shaded region). In the lower frequency range 
^ci(N) < uj < uj c i(1) (when do < d c ) and u c i(N) < uj < uj c2 (N) (when do > d c ), the modes 
are localized at the right hand side of the chain, and are called red plasmonic gradons (RG) 
(crossed shaded region). In the higher frequency range uj c2 (N) < uj < uj c2 (1) (when d < d c ) 
and uj c i(l) < uj < uo c2 {l) (when do > d c ), the modes are localized at the left hand side of the 
chain, and are called blue plasmonic gradons (BG) (tilted dotted line shaded region). For 
interparticle spacing do > d c , the modes become localized in the middle parts of the chain 
and covers a frequency band uj c2 (N) < uj < uj c i{1) (wavy shaded region). We call these 
modes blue-red plasmonic gradons (BRG) which are of great significance for the oscillations 
described below. Since the evolution of incoming wave depends on the initial conditions, we 
need to do the following wavepacket dynamic analysis. 

IV. WAVEPACKET DYNAMIC ANALYSIS OF VARIOUS OSCILLATIONS OF 
COUPLED PLASMON WAVES 

Once the complete set of gradon eigenmodes |0) has been obtained, we can construct an 
initial state to perform time domain simulations on the "wavepacket" dynamics. Let us take 
an initial wavefunction, 
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where ko is the vacuum wavenumber determining the central frequency and the direction of 
propagation of the wavepacket. The intensity profile \ip{x, 0)| 2 has a Gaussian distribution 
centered at x with spatial width a x . We then expand the initial wavefunction in terms of 

10), 

|^(O)) = ^A n ,|0 n ), (5) 

n 

where A n = (</> n |^(0)) is the constituent component of the initial wavepacket. The intensities 
\A n \ 2 of the various initial wavepackets have been shown as insets in Fig. 2. The peak of 
\A n \ 2 corresponds to the contribution from a dominant component of frequency u n , which 
is in agreement with the central frequency obtained from the dispersion relation u(k) at ko. 
Thus the subsequent wavefunction at time t is, 

|V>(t)} = 5>„|0 n } e ^. (6) 

n 

When the size N becomes large, one can replace the sum by an integral over the spectral 
components in Eqs. (5) and (6). The evolution of the wavepacket intensity \ip(x, t)\ 2 can 
be used to illustrate the various oscillations, including BO, breathing- wave-like (BW) os- 
cillations and other kinds of motions. We will show that combination of different localized 
modes (gradons) lead to different dynamic evolutions of the wavepacket. 

Here a few cases will be addressed to illustrate the correspondence between the various 
gradons and the oscillation of wavepackets. The occurrence of various oscillations depends on 
the localization extents of modes forming the initial wavepackets. Since the phase diagram 
as sketched in Fig. 1 explicitly reflects the spatial extension of the various gradon modes 
that the graded plasmonic chain can sustain, it can help us to have an instant judgment that 
whether an initial wavepacket contains a component that can reach either the left or the right 
end of the chain, where reflection occurs. Such kind of wavepackets shall not undergo Bloch 
oscillation or breathing-wave-like oscillation. A necessary condition for occurrence of BO or 
BW oscillation is that the initial wavepacket only consists of blue-red gradons, because BRG 
are localized in the middle part of the chain. We have marked the constituent components 
(intensity versus frequency) of the different initial wavepackets at certain do on the phase 
diagram (see Fig. 1). The corresponding dominant gradon mode of the initial wavepackets 
is shown in the insets, i.e. plots of |^(x)| 2 versus x. As expectations, in the BRG region, 
BO can occur. Depending on the initial wavepacket, breathing-wave-like (BW) oscillations 
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can also occur. The differences are that the initial wavepacket that undergoes BO has a 
larger spatial width (narrow bandwidth) while the initial wavepacket that undergoes BW 
oscillation has a relatively smaller spatial width (broad bandwidth). In the RG region, the 
initial wavepackets consisting of red gradons must be reflected at the right (red) end of the 
chain, such kind of dynamic evolution is denoted by RR (right end reflection). While in the 
BG region, left (blue) end reflection denoted by LR occurs. Finally, for the U region, an 
initially extended wavepacket can be reflected by both ends of the chain, which is denoted 
by LRR (left and right ends reflection). 

To show more clearly the above analysis regarding gradon localization and the occurrence 
of BO, we compose the contour plots of the intensity profile \ip(x, t)\ 2 in the position-time 
(i.e., x-t) domain for various initial wavepackets. The white color and red color indicate 
the very strong and relatively strong intensity, respectively, and the black color means the 
intensity is weak or zero. Figure 2(a) shows the evolution of \if)(x, t)\ 2 for d = 5.5a, a x = 5, 
and ko = 0.87T /do- The wavepacket exhibits an oscillatory motion: the mean position shows 
a periodic time-dependence while the width is nearly constant. This is a typical plasmonic 
BO process. More interestingly, in Fig. 2(b), we show the case of do = 5.0a, a x = 0.2, and 
ko = 0.87r/do. Now the wavepacket's width shows a periodic time-dependence but the mean 
position is nearly fixed at the initial place. This is a breathing-wave-like (BW) oscillation. 
Both initial wavepackets for BO and BW oscillation are only formed by blue-red gradons 
which are normal modes localized in the middle part of the graded chain. The typical modes 
at central frequencies uj = 0.243cj p (BO) and uj = 0.245u>p (BW) are shown in the insets of 
Fig. 1. However, the spatial widths of the two initial wavepackets are significantly different. 
For BO, the spatial width is larger than that of breathing- wave-like oscillation: the spatial 
width is a x = 5 for BO and a x = 0.2 for BW. 

We can estimate the period of Bloch oscillation. In semi-classical theory, 8 the equation 
of motion of the Bloch wavenumber k reads k = —du/dx = f, where / is almost constant 
for BO. The time taken for k to change by 2n/do is defined as the period of BO, denoted as 
T ■ 

2vr 2nN(2e l + c + Q 3 / 2 

Here the simplification is obtained by the fact that the average of cos k vanishes in a period. 
The plasmon resonant frequency 23 is u p = 1.72 x 10 16 rad/s. Thus the period of BO in 
Fig. 2(a) is about 11.8 ps. This is in agreement with the previous work for optical BO. 8 For 
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non-Bloch oscillations, the time between two neighboring reflections is defined as the period 
T R . 

Despite of the oscillation dynamics, we have also investigated the non-Bloch oscillations, 
e.g. wave reflection from the left end, the right end, or from both ends. Figure 2(c) and 
(d) show the corresponding situations for reflection from the left and right end, respectively. 
In Fig. 2(c) (do = 3.5a, a x = 5, ko — 0.97r/do), the initial wavepacket is constructed by 
linear combination of blue gradons only, which are localized modes residing at the left 
hand side of the chain (as shown in the inset named as LR in Fig. 1), this wavepacket can 
be reflected by the left end of the chain and can not reach the right end. In Fig. 2(d) 
(d = 4.0a, <j x = 5, ko = 0), the initial wavepacket is only formed by red gradons which are 
modes confined at the right hand side of the chain (as shown in the inset named as RR 
in Fig. 1), the wavepacket is reflected at the right end of the chain and can not reach the 
left end. If the components of initial wavepacket all fall into the extended modes region 
(see Fig. 1), the wavepacket can reach both ends and be reflected in multiple fashion before 
eventually spreading across the whole chain (the contour plot is not shown here). 

The different features of various oscillations can also be demonstrated in the contour 
plots of \ip(x, t)\ 2 in reciprocal position-time (i.e., k-t) domain. For BO or BW, k varies 
periodically in the range [— ir, ir] , while for other kinds of oscillations, k varies only in part 
of the range [— ir, n]. There are two reasons why we only show the contour plots of \ip(x, t)\ 2 
in real spatial position-time (i.e., x-t) domain. One is the contour plots in real space can 
show the obvious features of BO, BW, or other kinds of oscillations, but those in k space 
only show the periodic variation of k with time, which is sensitive to the initial value of k. 
The second reason is that the difference between BO and BW is more obvious in real space 
than that in k space. 

V. EVOLUTION OF WAVEPACKETS IN DAMPING CASE 

To be more realistic, we shall take into account the loss in the metal nanoparticles. 
For the damping case (r ^ 0), the equation of motion is modified by simply replacing 
uj 2 with uj{uj + ir) that is used in Eq. (3). Thus the eigenvalues become complex- valued 
Cj = u 0n + i-f n , where 7 n = T/2 and u 0n = to n (l — T 2 /Su 2 ). Thus the wavepacket is 
damped by an overall factor e~ vt l 2 . Figure 3 (d = 5.5a, a x = 5, k = 0.87r/d ) shows 
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the evolution of wavepackets for undamped (T = 0, solid lines) BO in an ideal chain and 
damped (T = 0.08u;p, dashed lines) BO in silver nanoparticle chains. At t = to (black lines) 
and t = t + 0.5Tb (blue lines), the wavepacket reaches at the right-most and left-most of 
the BO region. The damped wavepackets almost overlap the undamp ones. We can observe 
that the wavepacket at the left-most is higher and narrower than that at the right-most. 
This is due to the width of wavepackets also oscillates with time, which is similar to the 
photonic BO in electrically modulated photonic crystals. 26 As time goes on, at t = to + 2.0Tb 
(red lines) and t = to + 4.0Tb (green lines), the intensity of damped BO becomes smaller 
than that of the corresponding undamp one, with the damping rate T/2. At t = to + 2.5Tb 
(orange lines) and t = £o + 4.5Tb (light green lines), the wavepacket of damped BO is higher 
and narrower than that of corresponding undamp one. The amplitude of oscillation (i.e., 
distance between the left-most peak and right-most peak) is smaller in damped case than 
that in undamp case. 



VI. DISCUSSION AND FURTHER WORK 



In this work, we concentrate on studying the correspondence of gradon confinements and 
the Bloch oscillation as well as non-Bloch oscillations. This is realized by applying a linear 
gradient potential in the periodic plasmonic chain. The linear gradient potential is brought 
by the graded host permittivity. If we replace the linear gradient potential field by using a 
harmonic trap, we will realize the dipole oscillation. For dipole oscillation, the momentum 
oscillates periodically in time at the bottom of the trap, the average of velocity is periodic in 
time, and the phase difference between neighboring particles remain locked. However, dipole 
oscillation is unstable in optical waveguides with transverse confinement. Long living dipole 
oscillations can be realized in photonic crystals with Kerr dielectrics by properly applying 
graded pump fields. 27 Transition between Bloch oscillations and dipole oscillations will be 
studied in further work. Besides, based on our understanding of transitions among different 
gradon modes and the correspondence between gradon confinements and various oscillations, 
we can realize the steering among various oscillations by tuning the frequency or the graded 
parameters. 



9 



VII. CONCLUSION 



In conclusion, we elaborated the localized coupled plasmon modes and studied the dy- 
namics of plasmonic wave in GPCs. We found that there exist obvious correspondences 
between gradon localization and occurrence of various oscillations. The condition for occur- 
rence of Bloch oscillation or breathing-wave-like oscillation is that the initial wavepacket is 
formed only by blue-red plasmonic gradons. If the initial wavepacket consists of other kinds 
of plasmonic gradon modes, e. g. blue (red) gradons, the wavepacket will be reflected by 
the left (right) end of the chain. The predicted correspondence between gradon localization 
of initial wavepacket components and the type of subsequent dynamic oscillations opens the 
possibility to design various optical devices. We employed the wavepacket dynamics here and 
exact numerical results were obtained. When the chain size becomes large, some alternative 
methods, e.g., semi-classical theory and finite-difference time-domain methods can be used. 
The evolution of damped Bloch oscillation in a realistic silver chain shows similar results as 
those of no damping case, except that the intensity gradually decreases. The quasi-static 
dipole interaction has offered a good approximation, which gives similar results as that of 
fully retarded interaction. 21 The marked phase diagram provides an effective tool for the 
experimental realization of Bloch oscillations and for structure design of graded plasmonic 
arrays. A clear understanding of the mechanism of Bloch oscillations and other kinds of 
motions enables us to make full use of the graded plasmonic arrays, as promising candidates 
for subwavelength optical circuit and optical storage devices. 
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Figure Captions 

FIG. 1: (Color online) Phase diagram for the graded plasmonic chain (N = 100, c = 1.0, e/ = 3.0, 
£ = 5.45, and r] = 0.73). There are four regions of different gradon modes in the phase diagram 
as follows, BG: blue gradons region, RG: red gradons region, BRG: blue-red gradons region, and 
U: unbounded modes region. The corresponding types dynamic motion types of wavepackets 
consisting of different gradon modes are named as follows, LR: left end reflection, RR: right end 
reflection, BW: breathing-wave-like oscillations, and BO: Bloch oscillation. The typical mode 
patterns of these wavepackets are shown in the insets. The extended mode in U region extends 
the whole system and is not shown here. Bloch oscillation and breathing wave-like oscillation 
correspond to the same type of mode in the blue-red gradon phase, but occur with different initial 
wavepackets. 

FIG. 2: (Color online) Dynamics of various wavepackets shown by the contour plots of \ip(x, t)\ 2 
on the x-t domain, (a) Bloch oscillation (BO) (do = 5.5a, a x = 5, ko = 0.87r/do, an d Tg = 11.8 
ps), (b) breathing-wave-like oscillation (do = 5.0a, a x = 0.2, ko = 0.87r/cio, and Tg = 11.8 ps), 
and (c) reflection from the left (blue) end (do = 3.5a, a x = 5,/co = 0.97r/do, and Tr = 7.3 ps) (d) 
reflection from the right (red) end (do = 4.0a, a x = 5, ko = 0, and Tr = 7.3 ps). 

FIG. 3: (Color online) Comparison of the evolution of wavepackets of a undamped BO (r = 0, 
solid lines) in the ideal nanoparticle chian and of a damped BO (r = 0.08u;j,, dashed lines) in the 
silver nanoparticle chains (do = 5.5a,a x = 5, ko = OSn/do)- 
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